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asserted identity of a candidate-person. In an enrollment 
phase, a reference PQRST heart action graph is provided or 
constructed from information obtained from a plurality of 
graphs that resemble each other for a known reference person, 
using a first graph comparison metric. In a verification phase, 
a candidate-person asserts his/her identity and presents a 
plurality of his/her heart cycle graphs. If a sufficient number 
of the candidate-person’s measured graphs resemble each 
other, a representative composite graph is constructed from 
the candidate-person’s graphs and is compared with a com- 
posite reference graph, for the person whose identity is 
asserted, using a second graph comparison metric. When the 
second metric value lies in a selected range, the candidate- 
person’s assertion of identity is accepted. 
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FIG. 8B To step 103 








U.S. Patent 


Dec. 30, 2014 


Sheet 7 of 10 


US 8,924,736 B1 



FIG. 8C 









U.S. Patent 


Dec. 30, 2014 


Sheet 10 of 10 


US 8,924,736 B1 



FIG. 11 







US 8,924,736 B1 


1 

BIOMETRIC SUBJECT VERIFICATION 
BASED ON ELECTROCARDIOGRAPHIC 
SIGNALS 

ORIGIN OF THE INVENTION 

This invention described herein was made in the perfor- 
mance of work under a NASA contract and by an employee of 
the United States Government and is subject to the provisions 
of Public Law 96-517 (35 U.S.C. §202) and may be manu- 
factured and used by or for the Government for governmental 
purposes without the payment of any royalties thereon or 
therefore. 

FIELD OF THE INVENTION 

This invention relates to verification of an asserted identity 
of a candidate person, using biometric indicia. 

BACKGROUND OF THE INVENTION 

One or more biometric indicia, such as fingerprints, voice- 
prints, retinal scans and facial features, are often proposed to 
be used to identify, or to authenticate the asserted identity of 
a user who seeks access to a given resource. Approximately a 
dozen different biometric indicia have been proposed, but 
implementation methods for some of these approaches have 
not been disclosed. Many of these biometric indicia are asso- 
ciated with inherent physiological characteristics of the 
user’s body. Another set of such indicia relate to what may be 
characterized as behavioral characteristics that partly reflect a 
learning or behavioral process and do not rely exclusively on 
purely physiological features. Use of one or more of these 
behavioral characteristics as a biometric indicium has 
received relatively little attention, in part because of the per- 
ceived difficulty of implementing a procedure to measure 
such a characteristic. 

An example of a new physiological biometric feature is a 
sequence of bioelectric signals associated with cycles of the 
heart. 

Cardiac muscle is myogenic and is capable of generating 
an action potential and depolarizing and repolarizing signals 
from within the muscle itself. An intrinsic conduction system 
(ICS), a group of specialized cardiac cells, passes an electrical 
signal throughout the heart, schematically illustrated in FIG. 
1 . The ICS includes a sino-atrial (SA) node, an atrio-ventrical 
(AV) node, the bundle of His, right and left bundle branches, 
and the Purkinje fiber. These components spread the depolar- 
ization waves from the top (atria) of the heart down through 
the ventricles. The autonomic nervous system modulates the 
rhythm, rate and strength of cardiac contraction. When the 
cardiac muscle fibers contract, the volumes within the two 
atrial or two ventricle chambers are reduced and blood pres- 
sure increases. The (smaller) atrial chambers receive blood 
from the veins and pump the blood into the (larger) ventricle 
chambers, which pump blood out into the major arteries. The 
heart cycle normally begins in the right atrial chamber, and 
spreads to the left atrial chamber and then to the two ven- 
tricles. The atrial contraction is followed by the ventricular 
contraction in each cycle. 

Simultaneous contraction of the large number of fibers in 
the ICS generates an electrical field that can be measured at 
the body surface using an electrocardiograph (ECG). This 
electrical signal includes a sequence of PQRST complexes, 
one of which is schematically illustrated in FIG. 2. Most 
PQRST sequences are not uniform. The time interval 
between two consecutive R signal peaks, referred to as an R-R 
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interval, corresponds to a heart pulse, with a rate that nor- 
mally lies in a range of 60-90 beats per minute (bpm). The P 
signal corresponds to atrial depolarization (right side, fol- 
lowed by left side); the larger QRS complex corresponds to 
5 depolarization of the ventricles and repolarization of the atria; 
and the T signal corresponds to repolarization of the ven- 
tricles. A weaker U signal occasionally appears on the chart, 
representing remnants of ventricular repolarization, but is not 
shown in FIG. 2. 

to A “wave” comprises a curve covering at least one complete 
component (P, Q, R, S and/or T). According to naming con- 
ventions accepted by most cardiology workers, a time incre- 
ment with a straight line amplitude extending between two 
consecutive signals, for example, from the end of an S wave 
15 to the beginning of an immediately following T wave, is 
referred to as a “segment.” A time increment that includes at 
least one wave, with a graph that is at least partly curved, for 
example, from the beginning of a Q wave to the end of an S 
wave, is referred to as an “interval.” Herein, a “wave,” such as 
20 a P wave, will refer to the curvilinear graph (only) portion of 
a time interval, not including the associated time segment. A 
“wave”, of a “segment,” and of an “interval” are illustrated in 
FIG. 3. 

A P-Q time interval, normally of length At(p-q)»l 20-200 
25 msec, represents conduction time from initiation of atrial 
repolarization until initiation of ventricular depolarization, 
which is conventionally measured from the start of the P wave 
to the start of the Q wave. 

Where the ICS is diseased or is affected by presence of 
30 Digitalis, the P-Q time interval may lengthen as the pulse rate 
decreases; a prolonged P-Q interval, substantially beyond 
200 msec in length, is often evidence of atrio-ventricular 
block. .An abnormally short P-Q interval, substantially below 
120 msec in length, is often associated with hypertension 
35 and/or with paroxysms of tachycardia. The P-Q interval can 
also be shortened where the impulse originates in the AV 
node, or other atrial locations, rather than in the SA node. 

The QRS time interval, normally of temporal length At(q-t) 
z, 50- 1 00 msec, represents conduction time from initiation of 
40 ventricular depolarization until the end of ventricular depo- 
larization. and includes spread of the electrical impulse 
through the ventricular muscle. The P wave signal is normally 
gently rounded, not pointed or notched, and has a temporal 
length«50-l 10 msec. One or more of the P, Q, R, S and/or T 
45 peaks may be inverted in FIG. 2, depending upon electrode 
placement. A QRS interval greater than about 120 msec in 
temporal length often indicates ventricular arrhythmia or a 
block of one of the bundles. 

Normally, an S-T segment amplitude is approximately 
50 equal to a P-Q segment amplitude. The amplitude of the S-T 
segment, relative to the baseline (e.g., elevated or depressed), 
and the shape of the T signal are often of interest. The T signal 
is normally rounded and slightly asymmetrical. Presence of a 
sharply pointed or grossly notched T signal may indicate 
55 presence of myocardial infarction (pointed segment) or of 
pericarditis (notched segment). 

In some subjects, a beat (a single PQRST cycle) is some- 
times missed. This arises from the particular physiology of 
that subject and has not (yet) been shown to arise unambigu- 
60 ously from the presence of high stress in that subject. 

In normal ECG practice, ten or more electrodes including 
a ground electrode, are attached to the subject at selected, 
spaced apart locations. A chart of each PQRST cycle is 
printed separately on a 1 mmxl mm grid, with 1 mm (hori- 
65 zontal) representing 40 msec (0.04 sec time increment) and 1 
mm (vertical) representing 0.1 millivolts (mV amplitude). An 
upper limit on amplitude is usually 20-30 nun (2-3 mV). 
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ECG analysis is generally limited to medical diagnostics 
and to comparison of shifts with the passage of time of ECG 
parameters of interest. No substantial work has been done 
applying the ECG results to other areas of interest, such as 
authentication of an asserted identity of a candidate-person, 
through analysis of selected ECG results to provide one or 
more physiologically based biometric indicia associated with 
the candidate-person. Nor has any substantial use been made 
of evidence of a malady such as myocardial infarction or 
pericarditis as a characteristic for verifying the identity of a 
candidate-person. 

What is needed is a method and associated system for 
measuring one or more (preferably several) statistical param- 
eters associated with PQRST cycle for a candidate person and 
authenticating, or declining to authenticate, the candidate- 
person’s asserted identify with a reference person, using a 
comparison of the measured statistical parameter values (bio- 
metric indicia) with corresponding reference parameter val- 
ues. Optionally, the comparisons should be cumulative so that 
the biometric indicia test can be made progressively more 
demanding, to minimize likelihood of commission of a type I 
error (positive result is false) and/or to balance the likelihoods 
of commission of a type I error and commission of a type II 
error (negative result is false) in these comparisons. These 
comparisons should also allow for possible changes with 
passage of time of PQRST cycle characteristics for a candi- 
date-person. Preferably, evidence of presence of a malady in 
a reference person should be available for biometric use in 
comparison of a candidate person with this reference person. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
method and associated system for authenticating, or declin- 
ing to authenticate, an identity asserted by a candidate-person 
(c-p). The PQRST cycles for two distinct persons are likely 
sufficiently different that a comparison will allow discrimi- 
nation between these persons, based on differences in models 
constructed for each of a reference person and a candidate- 
person. Further, if a reference person is known to have a 
cardiac-related anomaly, absence of this anomaly in the mea- 
sured PQRST complex for the candidate-person strongly 
indicates that the candidate-person is not the reference per- 
son, and conversely. Presence of this anomaly in the reference 
person and in the candidate-person supports, but does not 
necessarily require, a conclusion that the candidate-person is 
the reference person. 

The original heart PQRST heart signal may include a 
superimposed non-zero base line signal f(t,base) that arises 
from voltage drift, noise and other undesirable effects. The 
original heart signals (reference and observed) are initially 
processed in order to “normalize” this signal by suppressing 
low frequency and high frequency components that are irrel- 
evant to the real heart signal. In one approach, low frequency 
signals, fisf(low) (e.g., f(low)=2 Hz), and high frequency 
components, f>f(high) (e.g., f(high)=42 Hz) are suppressed 
or removed using Butterworth or other suitable digital filters 
of tenth order. FIG. 4 schematically illustrates ideal fre- 
quency filtering according this approach. This processing 
removes slow drift and high frequency noise and provides a 
“normalized” heart signal for which the observed and refer- 
ence heart signals are reproducible and can be reliably com- 
pared with each other. FIGS. 5 and 6 illustrate PQRST time 
graphs of an original unfiltered heart signal and a filtered 
signal with low frequency components and/or high frequency 
components removed. 
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With removal of an undulating base line signal, the result- 
ing filtered signal in FIG. 6 appears to follow an approxi- 
mately constant zero base line, as desired. One can achieve 
somewhat the same results by forming an average baseline 
5 signal, using a moving time window and subtracting this 
average baseline, whose amplitude will change slowly with 
passage of time. One advantage of use of a digital filtering 
approach is that removal of the base line contributions is 
automatic and occurs as each succeeding heart signal is 
to received. Base line contributions are preferably removed 
before the procedures in FIGS. 7A/7B and 8A/8B/8C for the 
reference and candidate-person graphs are implemented. 

In a first embodiment (enrollment phase), the system 
receives a sequence of heart electrical action cycles for a 
15 known reference person. These cycles are segmented and 
recorded, and a set SI of N graphs, numbered n=l, . . . , N, of 
these reference cycles is formed. K such graphs are randomly 
chosen (3<K<N), numbered k=l, . . . , K. A pair of graphs, 
v(t;k=kl ) and v(t;k=k2), is identified for which a first graph 
20 comparison metric p 1 (kl ;k2) is a maximum for graphs drawn 
from the chosen set, and a composite graph, denoted v(t;kl 2), 
is formed using data from the graphs v(t;kl) and v(t;k2). 

A set S2 of graphs (n) among the N graphs is identified for 
which a metric value pl(kl2;n) is at least equal to a selected 
25 threshold value (thrl). The threshold value (tlxrl) may be 
chosen so that the set S2 includes a selected fraction g of the 
N graphs that resemble each other, with g preferably lying in 
a range of about 0.3sg<0.99; a higher choice of the threshold 
value (tlxrl ) will lead to a correspondingly lower value of the 
30 fraction g. From each graph in the set S2, a feature vector FV, 
discussed in Appendix C, is constructed, representing an 
amplitude and corresponding time position of local maxima 
and minima of the graph. From the resulting set of feature 
vectors, a Gaussian mixture model, denoted GMM(S2;ref), is 
35 constructed and used as a reference heart cycle model, repre- 
senting the known reference person, relying upon the fraction 
g of graphs that resemble each other. Preferably, a Gaussian 
mixture model, constructed in accord with the article entitled 
“Maximum Likelihood from Incomplete Data via the EM 
40 Algorithm,” by Dempster, Laird and Rubin (Trans. Royal 
Statistical Soc., Vol. 39, Series B, pp. 1-38 (1977), incorpo- 
rated by reference herein, is used here. 

In a second embodiment (verification phase), the system 
subsequently receives a sequence of heart cycles and associ- 
45 ated graphs for a candidate-person (c-p), who asserts an iden- 
tity with a reference person that is to be verified or refuted. A 
sequence of N' graphs v'(t;n';c-p), numbered n'=l, . . . , N', is 
presented by the c-p, and a set S 1 ' of K' such graphs v'(t ;k' ;c-p ) 
(k’=l , . . . , K';3<K'<N') is randomly chosen. A pair of graphs, 
50 v'(t;k'=kl';c-p) and v'(t;k'=k2';c-p), is identified for which a 
metric value p 1 (kl ’;k2’) is a maximum for graphs drawn from 
the chosen set 51 ; and a composite graph, denoted v'(t;kl2'; 
CP), is formed using data from the graphs v'(t;kl') and v'(t; 
k2'). 

55 A set S2' of graphs (n') among the N' graphs is identified for 

which a metric value pl(kl2';n') is at least equal to a selected 
threshold value (thr2). The threshold value (thr2) may be 
chosen so that the set S2' includes at least a selected fraction 
g' of the N' graphs that resemble each other, with g' preferably 
60 lying in a range 0.3<g'<0.99. From each graph in the set S2, a 
feature vector FV 1 is constructed representing an amplitude 
and corresponding time position of local maxima and minim a 
of the graph. From the resulting set of feature vectors, a 
Gaussian mixture model, denoted GMM(S2';rep), is con- 
65 structed and used to represent the c-p. A second metric, p2, 
referred to as a GMM comparison metric, is provided and 
applied to compute a second metric value, p2{GMM(S2;ref); 
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GMM(S2';c-p)}, which compares the reference GMM model 
with the c-p GMM model. The candidate-person’s asserted 
identity is accepted as likely to be valid if (and, optionally, 
only if) the second metric value p2{GMM(S2;ref);GMM 
(S2';c-p)} is at least equal to a third selected third threshold 5 
value (thr3). If the second metric value p2{GMM(S2;ref); 
GMM(S2';c-p)} is less than a fourth threshold value (thr4), 
with (thr4)<(thr3), acceptance of the candidate-person’s 
asserted identity is optionally declined; or one or more further 
tests of the candidate-person’s identity are applied. 10 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates a human heart. 

FIG. 2 schematically illustrates a representative PQRST 15 
electrical signal complex generated in association with a 
human heart. 

FIG. 3 illustrates “waves”, “segments” and “intervals.” 

FIG. 4 illustrates ideal frequency filtering of an original 
heart signal 20 

FIGS. 5 and 6 graphically illustrate a result of removal of a 
base line contribution to an original, unfiltered heart signal. 

FIGS. 7A/7B and 8A/8B/8C are flow charts illustrating an 
enrollment phase and an identity verification phase for prac- 
ticing the invention. 25 

FIG. 9 illustrates estimation of an anchor central point 
(acp). 

FIG. 10 illustrates contraction and expansion of time scale 
with heart rate. 

FIG. 1 1 is a flow chart of a procedure to construct a feature 30 
vector. 

DESCRIPTION OF THE INVENTION 

FIGS. 7A/7B are a flow chart illustrating an embodiment of 35 
an enrollment phase for practicing the invention. In step 70, a 
counting index c is initially set equal to 0 (c(initial)=0). In step 
71, the system provides or receives and stores a sequence of 
heart electrical action cycles (sometimes referred to as 
“PQRST signatures”) to create a model for a reference person 40 
whose identity is known. Preferably, each cycle covers only 
one PQRST heart electrical action time interval, and the 
cycles do not overlap, except optionally at the cycle end 
points. In step 72, the system forms a time dependent graph 
v(t;n) for each of N (N>3) recorded heart electrical action 45 
cycles, numberedn=l , . . . , N. In step 73, the system randomly 
chooses a set of K of the heart cycle graphs (3<K<N), num- 
bered k=l , 2, . . . , K, and compares each chosen graph v(t;k) 
with each other chosen graph, using a first selected, numeri- 
cal-valued graph comparison metric pi. 50 

In step 74 the system identifies a pair, v(t;kl ) and v(t;k2), of 
the chosen reference graphs that has a metric value pi (kl,k2) 
that is as larger as or larger than any other metric value 
pl(k',k") among the K chosen reference graphs. In step 75, 
the system forms a first composite heart cycle graphv(t;kl2), 55 
using data from the graphs v(t;kl ) and v(t;k2). 

In step 76, metric values pl(kl2;n) are computed for all 
graphs v(t;n) (n=l, . . . , N) in the set of N graphs. In step 77, 
the system queries whether at least a selected fraction g 
(0.5<g<l) of these metric values satisfy pl(kl2;n)>(tlirl), 60 
where (thrl) is a selected first threshold value. Where the 
answer to the query in step 77 is “yes,” a subset S2 of all 
graphs v(t;nl) (nl=l, 2, 3, ... ) in the set SI is identified for 
which pl(kl2;nl)>(tlirl), in step 78. Choice of a higher 
threshold value (thrl) will lead to a correspondingly lower 65 
value of the fraction g. Optionally, the threshold value (thrl ) 
may be chosen in a range of about (thrl )=(0.3-0.99) pl(kl, 
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k2). Alternatively, the criterion for graph resemblance may be 
pl'(kl2;nl)<(tlirl'), for another selected threshold value 
(thrl 1 ), with pT(kl,k2)<(thrl'), for an alternative first graph 
comparison metric pi'. 

For each graph in the set S2, a feature vector FV is con- 
structed, for example, as in Appendix C, representing the 
amplitude and corresponding time position of local maxima 
and minima of the graph v(t;nl ) in the set S2, in step 79. From 
the resulting set of feature vectors, a sum of Gaussian mixture 
model (Appendix B), denoted GMM(S2;ref), is constructed 
and optionally stored, in step 80, and used in step 81 as a 
reference heart cycle model, representing the known refer- 
ence person, incorporating data from the graphs v(t;nl) inS2. 
Preferably, a Gaussian mixture model, constructed in accord 
with the article entitled “Maximum Likelihood from Incom- 
plete Data via the EM Algorithm,” by Dempster, Laird and 
Rubin (op. cit.), is used here. 

Where the answer to the query in step 77 is “no,” the 
system, in step 82, increments the counting index c (c— »c+l). 
In step 83, the system inquires if the incremented index c 
satisfies c>c(max), where c(max) is a selected maximum 
index value (e.g., c(max) between 20 and 100). If the answer 
to the query in step 83 is “no,” the system returns to step 73 
and repeats steps 73-83. If the answer to the query in step 83 
is “yes,” the system, moves to step 84 and interprets this 
condition as indicating that it is impossible to create a refer- 
ence model for this reference person. 

Appendix B sets forth a procedure for construction of a 
Gaussian mixture model, consisting of a sum of individual 
Gaussian functions, each having a weight or amplitude, a 
temporal mean or location of an extreme value (maximum or 
minimum), and a characteristic width relative to the peak 
value (proportional to standard deviation of a Gaussian func- 
tion). A Gaussian mixture is a sum of M Gaussian functions 
so that the total (maximum) number of parameters or “fea- 
tures” is 3M. A feature vector is formed from these features 
for each of the graphs v(t;nl), and a model, relying on data 
from the feature vectors, is formed. A smaller number of 
features (<3M) may be used to form a feature vector. One 
result is a probability density function p(xl3M) for feature 
vector components that correspond to amplitude values, 
mean values and standard deviation values for M Gaussian 
functions that together represent the heart action graphs in 
different time sub-intervals of a heart action cycle. 

FIGS. 8A/8B/8C are flow charts illustrating an embodi- 
ment of an identity verification phase for practicing the inven- 
tion. In step 90, the system sets an initial value of a counting 
index c' (e.g., c'(initial)=0). In step 91, the system receives an 
assertion of an identity for a candidate-person. In step 92, the 
system provides or receives and records a sequence of heart 
electrical action cycles or PQRST signatures for the candi- 
date-person. In step 93, the system forms a heart cycle graph 
v'(t;n') for each of a set of N' (N'>3) of the reference heart 
electrical action cycles, numbered n'=l , . . . , N', presented by 
the candidate-person. Preferably, the cycles do not overlap, 
except optionally at the cycle end points. 

In step 94, the system randomly chooses a set of K' of the 
c-p heart cycle graphs (3>K'>N'), numbered k'=l , 2, . . . , K', 
and compares each chosen graph v'(t;k') with each of the other 
chosen c-p graphs, using the first graph comparison metric 
pi. In step 95, the system identifies a pair of chosen c-p 
graphs, v'(t;kl') and v'(t;k2') among the chosen graphs for 
which the first metric value pl(kl';k2') is at least as laige as, 
or larger than, any other metric value pl(k'",k"") among the 
K' chosen graphs. In step 96, the system forms a composite 
c-p heart cycle graph v'(t;kl2'), using data from the c-graphs 
v'(t;kl') and v'(t;k2'). 
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In step 97 , metric values pl(kl2';n') are computed for all 
graphs v'(t;n') (n'=l, . . . ,N') in the set of N' c-p graphs. Instep 
98 , the system determines if at least a selected fraction g 
‘(0.5<f<l) of these metric values satisfy pl(kl2';n')>(thr2), 
where (thr2) is a selected second threshold value. Where the 
answer to the query in step 98 is “yes,” a subset S2' of all c-p 
graphs v'(t;nl') is identified for which pl(kl2';nl')>(thr2), in 
step 99 . Choice of a higher threshold value (thr2) will lead to 
a correspondingly lower value of the fraction g 1 . Optionally, 
the threshold value (thr2) may be chosen in a range of (thr2)= 
(0.3-0.99) pl(kl',k2'). Alternatively, the criterion for graph 
resemblance may be pr(kl2';nT)<(thr2'), for another 
selected threshold value (thr2'), with pl'(kl',k2 , )<(thr2'), for 
an alternative first graph comparison metric pi 1 . 

Where the answer to the query in step 98 is “yes,” the 
system, in step 100, creates or provides a representative c-p 
heart cycle graph, v'(t;rep)=v'(t;kl2'), ^ optionally stores 
the representative graph v'(t;rep) for this candidate-person. 

In step 101 , the system compares the representative graph 
v'(t;rep) for the candidate-person and a reference graph v(t; 
ref) associated with the known reference person whose iden- 
tity is asserted by the candidate-person, using a second, 
numerical-valued graph comparison metric p2(rep;ref). The 
first and second metrics, pi and p2, may be different or may 
be the same. Where p2(rep;ref)>a selected second metric 
threshold value (thr2), the system, in step 102, interprets this 
condition as indicating that the candidate-person is likely to 
be the reference person whose identity is asserted by the 
candidate-person. Where p2(rep;ref)<a selected third metric 
threshold value (tlir3), with 0<(thr3)s(thr2), the system inter- 
prets this condition, in step 103 , as indicating that the candi- 
date-person is not likely the reference person whose identity 
is asserted by the candidate-person. 

Where the answer to the query in step 98 is “no,” the 
system, in step 104 , increments the index c' (c'-»c'+l ). In step 
105 , the system inquires whether the incremented index c' 
satisfies c'>c'(max), where c'(max) is a selected maximum 
index value (e.g., c'(max) between 20 and 1 00). If the answer 
to the query in step 105 is “no,” the system returns to step 93 
and repeats steps 93 - 105 . If the answer to the query in step 
105 is “yes,” the system, in step 106 , interprets this condition 
as indicating that it is impossible to create a representative 
graph for this candidate-person. Optionally, the candidate- 
person’s assertion of identity may be rejected. 

The first graph comparison metric pi (or pi') is optionally 
configured to compare values at corresponding times for a 
chosen pair of heart cycle graphs, indexed as (k'",k"") (k'", 
k""=l , . . . , K; k'"=k""), of the K graphs, or of the graphs kl 2 
and n, at each of a sequence of spaced apart time points, t=t m , 
at which the corresponding pair (k'",k"") of heart cycle values 
v(t m ;k) were provided and/or received, measured from a 
selected initial time value for each of the pair (k"',k"") of heart 
cycles. A suitable first metric pi is a normalized correlation 
function, defined as 

pl(i:"',i"'0=corr(i"',i"")/{coir(*"';i"')corr(i"";i"'9 

} 1/2 ( 1 ) 


« (2) 
coir (k”' , k"") = V >v(f OT )v(r OT ;A:" , )’*v(f m .A:"") 5 


where * indicates a suitable conjugation operation and w(t„,) 
is a selected non-negative weight value. In this approach, one 
requires that 


8 

in order that the graphs corresponding to {v(t m ;k"')}„, and 
{v(t m ;k"')} m , sufficiently resemble each other. 

A suitable version of an alternative first metric pi' is a sum 
of differences, 

5 

M (4) 

pT(*"', *"") = Yj w'lOKf™, k"")l V(k"’) - v(r„ ; *"")/ V(*"")K . 

m= 1 

10 « (5) 

v(k) = Y M'™;*)!. 


where w'(t m ) is a selected non-negative weight value (e.g., 
15 uniform weights), q' is a positive exponent number, and 
{v(t m ..k)} m , is a selected sequence of heart cycle values, 
associated with a known reference person. In this approach, 
one requires 

20 pl'(i"'i"")-(thrl')=Z)2sO (6) 

in order that the graphs corresponding to {v(t„,;k")} m and 
{v(t m ;k")} m , sufficiently resemble each other. The two 
inequalities in Eqs. (3) and (6) can be expressed in a single 
inequality, p 1 (k';k")-(thrl )=D1 , with D1 >0 or D1 <0, respec- 
23 tively. 

A second metric value p2 may also be defined according to 
Eqs. (l)and (2), or alternatively according to Eqs. (4) and (5), 
with possibly-different choices of weight values w"(t m ) and/ 
or possibly a different choice of an exponent index q', 

30 

p2(rep;ref)=corr(ref;rep)/{corr(ref;ref)corr 

(rep;rep)}‘« (7) 


35 M (8) 

con irep\ ref) = Y ref ) * v(t m ; rep). 

m= 1 

Where 

40 

p2(rep;ref)-(thr2)=£>2s0 (9) 

is satisfied for a selected second metric threshold value (thr2), 
this condition is interpreted as indicating that the candidate- 
4 _ person is likely to be the reference person. 

A suitable version of an alternative second metric p2' is a 
sum 


M ^ (10) 

30 p2'{rep; ref) = ’Y w ’ ^ ref) / V(ref) - v'(l m /ep) j V'(rep )\ q ' , 

m= l 

M ( 11 ) 

V(ref) = Y\ v ('m-,ref)\, 

m= 1 

55 

u ( 12 ) 

V’(rep) = Y b'h™; rep ) |, 


60 where w’(t m ) is a selected non-negative weight value, q' is a 
selected positive number, and {v(t m .rep} m and {v(t„,.rep)} m 
are selected sequences of heart cycle values, the latter being 
associated with the candidate-person. Where 

^ p2 , (rep;ref)-(thr2 , )=D2'sO (13) 

is satisfied for a selected third metric threshold (thr3), with 
0<(thr3)<(thr2), the system interprets this condition as indi- 


pl {k"',k”")-(fki\)=D\ £0. 


(3) 
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eating that the candidate-person is not likely the reference 
person whose identity is asserted by the candidate-person. 

A reference heart cycle graph can be formed using data 
from a pair of graphs, (k=kl, k=k2), that have a maximum 
first metric value pl(kl,k2). One suitable reference graph is 5 


pi (kl,k2)=max k tj e f'{pl(k , ",k'"')} (14) 

v(r;ref)=Avg il-i2 {v(r;W),v(rA2)}, (15) 

where Avg is a point-by-point averaging operation, such as io 

Avg{v(f;Al),v(?;A2)}={v(r;Al)+v(f;A2)}/2, (16-1) 

Avg{v(f;£l),v(f;A2)}=sgn{v(r;A'l)+v(r;A2)}lv(r;A'l)-v(f; 

kl) I" 2 , (16-2) 

15 

Avg{v(t;kl),v([;k2)}=v([;kl )+v(t;k2)-sgn { v{t;k\ )+v(f; 

k2)}\v(t;k\)-v(t;k2)\ m . (16-3) 


The threshold value (tlirl) may be in a range, such as (0.3- 
0.99) Max{pl(kl,k2)}. 

Another suitable reference heart cycle graph can be formed 20 
from a pair of graphs, (k=k3,k=k4), that have a minimum first 
metric value pl'(k3,k4), with the same expression v(t;ref,) as 
in Eq. (15) and in the examples in Eqs. (16-1), (16-2) and 


pi \k3,kA)=mm k > ^»{pl \k',k ,r )} 

(17) 

v(f;ref) = Avg { v(t;k3) , v(t;k4) } , 

(18) 


The threshold value (tlxr 1 ') may be a selected fraction (e.g., 

1-3 times the minimum value pl'(k3,k4). The considerations 30 
in Eqs. (14)-(18) may also be applied to formation of the 
representative graph values v(t;rep). 

The representative candidate-person heart cycle graph may 
be computed by averaging the heart cycle graphs (k"'=kl', 
k"'=k2') for the candidate-person. The heart rate-based cor- 35 
rections of Bazett/Fredericia (At(q-t)) and Danter-Carrothers 
(At(p-r)), discussed in Appendix A, are optionally applied to 
provide an HR-corrected average heart cycle graph for the 
known reference person and for the candidate-person, as 
defined on a time interval having a selected length (e.g., 40 
At=600 msec) within the graph. A biometric feature vector is 
extracted from the HR-corrected average heart cycle graph 
for the candidate-person, as discussed in Appendix C, of 
dimension d=18, or smaller if desired. 

A sequence of PQRST cycles associated with a reference 42 
person or a candidate-person will vary from one cycle to 
another, which requires a standardized determination of 
where a cycle “begins.” One approach focuses on a beginning 
of a P-wave, which is taken to be a time near the beginning of 
a cycle where the signal amplitude of the P-wave portion first 5 0 
diverges significantly from a zero-amplitude (isoelectric) 
baseline. One problem encountered here is that of distin- 
guishing between a non-zero P-wave amplitude that is merely 
noise, from a non-zero signal amplitude that is the begimiing 
of a true P-wave. 

One approach adopted here focuses, not on a “beginning” 
of a P-wave, but on an analog of a level crossing, sometimes 
referred to as a zero crossing. A time derivative of the PQRST 
signature signal amplitude U(t;k) is squared for each of a 
sequence of times {t„(k)}„ for each cycle m, and a time, 60 
t=t(k;max) is determined for each cycle for which U(t)) 2 is a 
maximum, 

t(kO;m£Lx)={t n (m)\(U(t) 2 is maximum at f=7(&0)}. (19) 

A squared level value L, such as 

L =/zmax{ I U(t) 2 }, (20) 


10 

(e.g., h=0.2 or 0.5 or 0.6) and a minimum first time value, 
t=Tl, and a maximum second time value, t=T2 (>T1), are 
estimated for which 


nzt(k0)*n, 

(21) 

{U{tf{t=T\)=U{i) 2 {t=Tl)=L, 

(22) 


as illustrated in FIG. 9; t=T 1 is the largest time value, less than 
t(k0), for which U(t) 2 (t=Tl) lias the value L, and t=T2 is the 
smallest time value greater than t(k0) for which U(t) 2 (t=T2) 
has the value L. The time value t=t(k0) will be located near, 
but not necessarily coincident with, a time value for which the 
R-wave achieves its largest amplitude. A mid-point time 
value, t=t(acp)=(Tl+T2)/2, will be referred to herein as an 
“anchor central point” (acp ) for a heart cycle and is optionally 
used in compensating for heart beat rates that differ from 
HR=60/min. 

A time interval including one heart cycle is now divided 
into two (or three) sub-intervals, relative to the acp for the 
cycle: 


7a.7(acp)-250 msec</s?(acp), 

(23A) 

76.7(acp)sf<r(acp)+350 msec. 

(23 B) 


Where three sub-intervals are provided for a cycle, the time 
division is often determined as 

71:f(acp)-250 msecsrsr(acp)-70 msec; 

72:?(acp)-69 msecsfsf(acp)+110 msec; 

73:?(acp)+l 1 1 msecs/s/(acp)+250 msec. 

The time t=t(acp) corresponds approximately to an aver- 
aged time for maximum slope of the R-wave component of a 
collection of cycles. The time interval la will include all of the 
P-wave (if present), all of the Q-wave (if present) and most of 
the R-wave that precedes t(acp), for most cycles; and the time 
interval lb will include the remainder of the R-wave, all of the 
S-wave (if present) and all of the T-wave (if present) for most 
cycles. The time intervals la and lb, 600 msec in total length, 
determine the time intervals for each heart cycle used in the 
computations of correlations, weighted differences and aver- 
ages fonned in Eqs. (1 )-(l 8). 

The representative heart cycle graph of values in each of 
the three time intervals, 11,12 and 13, for a PQRST signature 
can be constructed, using a Gaussian mixture (sum) of M=2 
or M=3 Gaussian functions, with each Gaussian, time depen- 
dent function having as many as three adjustable parameters: 
for example, maximum amplitude, time corresponding to 
mean value (maximum magnitude), and standard deviation. A 
procedure for determination of the parameters for the Gaus- 
sian mixture is set forth in Appendix C. 

A feature vector, having parameter values drawn from the 
adjustable parameters of the Gaussian mixture, will have Q 
vector components, with Q<18; in one instance, Q=14 is 
preferred. 

Characterization. 

The foregoing procedure includes: (i) an enrollment phase, 
wherein a reference graph) s) of heart electrical action for one 
or more known reference persons is created for subsequent 
comparison with a representative graph of heart electrical 
action presented by a candidate-person; and (ii) an identity 
verification phase, wherein the representative graph pre- 
sented by the candidate-person i s compared with one or more 
extant reference graphs, to confirm or refute the candidate- 
person’s assertion of identity with a reference person. The 
approach is strongly non-parametric; a reference graph is 
constructed using two or more graphs for the reference person 
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that resemble each other and that resemble a large fraction of 
other graphs associated with the known reference person. The 
system takes account of the possibility that no composite 
graph, formed from one or two graphs that most closely 
resemble each other, may adequately characterize heart elec- 
trical action for the known reference person. In the verifica- 
tion phase, heart electrical action graphs presented by the 
candidate-person are examined to estimate whether a com- 
posite graph, formed from the graphs presented, adequately 
characterizes heart electrical action for the candidate-person. 
Identity of the candidate-person with the reference person is 
confirmed if (1) a graph comparison metric, applied to the 
reference and representative graphs, has a value in a selected 
range and (2) the reference person graph and the candidate- 
person graph are both sufficiently well characterized in the 
respective populations from which they are drawn. Use of 
non-parametric comparisons (i) ensures that no heavy 
emphasis is imposed on having adequate (non-sparse) statis- 
tical data for such purpose, and (ii) facilitates comparisons 
where the graphs presented by the candidate person include a 
relatively large fraction of PQRST signatures with what 
appears to be pathological behavior. 

Appendix A 

Heart Cycle Correction Based on Overall Heart Rate 

In 1920, Bazett and, separately, Fredericia proposed an 
empirical relation between resting heart rate HR and a length 
of a QT interval. A corrected QT interval length, denoted 
AQTc, is observed to be approximated as 

AQT corr (msec)=Ag7tl{HR/60} 1/2 (Bazett). (A-l) 

AQT corr (msec)=A£)7D{HR/60} 1/3 (Fredericia). (A-2) 

where AQTO is the observed and uncorrected QT interval 
length for a subject with resting heart rate HR, normally in a 
range 40-120 beats per minute (bpm). The corrected QT 
length is reduced, relative to the uncorrected QT length, when 
the heart beat rate is less than 60 bpm, and inversely. FIG. 10 
illustrates a contraction in time scale, for HR=120 bpm (>60 
bpm) and illustrates an expansion in time scale, for HR=40 
bpm (<60 bpm), for the QT time interval. 

An analogous correction, APRc, expressed as a linear 
regression, for a PR length with a resting heart rate differing 
from 60 bpm has been proposed by Danter and Camithers 
(Brit. Jour. Clinical Pharmacol. Vol. 30(1 990) pp 490-492) as 

APR corr (msec)=l 95. 36-0.3 84-HR=A/ > i?0-{AF7t cor / 

APR0} (A-3) 

FIG. 10 also illustrates contraction of the PR time interval and 
expansion of the time interval, for HR=120 bpm and for 
HR=40 bpm. Where HR differs from 60 bpm, the time inter- 
val I and the time interval II are each expanded or contracted 
to corrected time intervals as shown, determined by the rela- 
tionships 


I (corr):-250-a (msec)s7s7(acp), 

(A-4) 

II (corr):7(cap)s7s350 p (msec)-, 

(A-5) 

a={l±PR corr J{iPm}, 

(A- 6) 

^{AQT^AQTO}, 

(A- 7) 


The (uncorrected) reference graph v(t;ref) and the (uncor- 
rected) representative graph v(t;rep) are replaced by corrected 
graphs with appropriately expanded or corrected time, 
according to 


12 

v(7;ref) co ^=v{r(acp)-a(7(acp)-7;ref)} (-250 (msec) 

srsr(acp)). (A- 8) 

v(7;ref) co ^=v{r(acp)+p(7-7(acp);ref)} (?(acp)sfs350 

(msec)), (A-9) 

v(7;rep) corr =v{7(acp)-a(?(acp)-7;rep)} (-250 (msec) 

ste?(acp)), (A10) 

v(7;rep) corr =v{ 7(acp ) + (5 (7- 7(acp );rep ) } (7 (acp) <X<. 3 5 0 
10 (msec)). (All) 

The corrected graph values v(t;ref) co , and v^rep)^,.,,, are to 
be used in place of v(t;ref) and v(t;rep), respectively, in the 
preceding. 

15 Appendix B 

Determination of Optimum Gaussian Mixture Model 
Parameters 

211 The representative heart cycle graph of values in each of 
the three time intervals, II, 12 and 13, for a PQRST signature 
can be constructed, using a Gaussian mixture (sum) of M 
Gaussian functions, with each Gaussian, time dependent 
function being represented as 

25 

G(f;a,;n 4 ,a)=a„[23ia„ 2 ] 1,2 exp{-(f-/j„) 2 /(2a„ 2 ), (B-l) 

where a m is a weight or amplitude (positive, negative or zero), 
|.i„, is a time value (mean value) selected for maximum mag- 
nitude for the heart electrical action cycle in the time interval 
3< 1 m=l , . . . , M, and ct„, is the associated standard deviation for 
the time interval m. For a given time interval i, the parameters 
a OT , p,„ and a m are chosen jointly to satisfy three minimization 
crit eria with respect to each of these three parameters of an 
error value 
35 



40 

where A ( (Q) is a measured amplitude of a variable of interest 
(e.g., heart electrical voltage developed) in the time interval 
m=l , .... M at a measurement time t=t A (k=l , . . . , K) and q 
45 is a selected positive number (e.g., q=l or 2 or /T7). The 
computed value of e is jointly minimized as 

with respect to each of the parameters a m , and a m , for 
50 example, by partial differentiation of e with respect to each of 
these parameters. The result is a set of 3N coupled homoge- 
neous equations to be solved simultaneously. 

A maximum likelihood (ML) method is preferably used to 
identify the parameters a m , \i m and a m associated with the 
55 individual Gaussian functions in Eq. (B-2). A useful discus- 
sion of application of an ML method is presented by Douglas 
Reynolds in “Speaker identification and verification using 
Gaussian mixture speaker models”, Speech Communication, 
vol. 17 (1995) pp 91-108. These computations are for each of 
60 the fraction g of graphs v(t;n) presented for the graphs (“re- 
semblance graphs”) that sufficiently resemble each other, 
according to the criterion pl(kl2;n)>(thrl ). Optionally, the 
fraction g of these resemblance graphs are divided into Q 
groups of such graphs (Q>2), and each group is subjected to 
65 a K-means clustering procedure. For further discussion here, 
it is assumed that Q=1 so that subsequent statistical process- 
ing is applied to one collection of the resemblance graphs. 
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After the optimum parameter values are determined for 
each of the resemblance graphs, each of the individual resem- 
blance graphs is approximated as 


v(r; n; approx ) = ^ a m [ \ 1/2 exp{-(r- p m ) 2 j (2cr^)} 


where M is the number of Gaussian functions used for each 
resemblance graph (number n). Each such Gaussian function 
has three parameters (weight or amplitude, temporal location 
of mean or peak amplitude, and standard deviation) so that 
3M parameters are used to approximate each resemblance 
graph (assumed to have a count of N). Each of the 3M param- 
eters corresponds to a feature. A 3Mx3M covariance matrix 
C=C(wl,w2) is calculated for each pair features (wl,w2) for 
all the Nf resemblance graphs, and a probability density func- 
tion 

p(x;3M)={2n) D/2 ||C|| 3/2 }- 1 exp{-(x-^)(Q- 1 (x-|i)' r /2}, (B-4) 

where x and p are a feature vector and a mean value vector, 
each having 3M dimensions. Occurrence of the inverse 
matrix C _1 in the exponent in Eq. (B-4) compensates for 
differences in dimensions and numerical values of the vector 
components x. 

Appendix C 

Creation of a Feature Vector for a Reference Person 

From each of the fraction g (or more) of the N graphs, one 
extracts a vector of features, of dimension D (to be specified). 
A Gaussian mixture model (GMM) is created, from a sum of 
M Gaussian functions (preferably, M=2 or 3), using an esti- 
mate-minimize algorithm corresponding to a known refer- 
ence person, with M weighting sets (a m , \i m , s„,) (m=l M) 

of parameters to be chosen to minimize an error function. 

FIG. 11 is a flow chart illustrating a procedure for con- 
structing a feature vector. In step 111, a time interval of length 
At(acp) (e.g., At(acp)=600 msec) is selected within each heart 
electrical action cycle for the majority of graphs for the 
reference person. Each such time interval comprises a first 
time interval II of length Atl (e.g., Atl = 180 msec), a second 
time interval 12 of length At2 (e.g., At2=l 80 msec), and a third 
time interval 13 of length At3 (e.g., At3=240 msec), where 
consecutive time sub-intervals do not overlap, except option- 
ally at an endpoint. 

In step 112 , a selected portion of each graph of each mea- 
sured heart electrical action cycle interval in the graphs is 
corrected (expanded or contracted), using the Bazett/Frederi- 
cia model for a At(q-t) section of each cycle and using the 
Danter-Carrothers model for a At(p-r) section of each cycle, 
based on observed, time-averaged heart rate (HR) for the 
reference person. 

In step 113 , the system then identifies a first local algebraic 
minimum amplitude, a first local algebraic maximum ampli- 
tude, and a second algebraic minimum amplitude, in that 
order, within each of the three sub-intervals (II, 12, 13), the 
time values for the corresponding minima, maxima and 
minima amplitudes, and the corresponding minima, maxima 
and minima amplitude values, respectively. This provides 1 8 
parameter values for each heart action cycle. 

In step 114 , this ordered sequence of 18 parameter values 
forms a feature vector FV of length d< 1 8 for each heart action 
cycle in the majority of graphs for the reference person. 
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Optionally, a feature vector FV (of length d<18) can be 
formed from the minimum and maximum amplitudes and 
corresponding time values discussed in the preceding, with- 
out application of the heart rate corrections for the Dt(q-t) 

5 and/or Dt(p-r) set forth in the Bazett, Fredericia and/or 
Danter-Carrothers articles. 

What is claimed is: 

1. A method of authenticating or declining to authenticate 

io an asserted identity of a candidate person (c-p), the method 

comprising: 

(i) initializing a count index c (c (initi al )=0) ; 

(ii) providing or receiving a sequence of N reference heart 
electrical action cycles, numbered n=l, . . . , N, of a 

15 known reference person, and forming a graph v(t;n) 
representing each provided heart cycle for the reference 
person; 

(iii) randomly choosing K heart cycle graphs, numbered 
k=l, 2, . . . , K (3<K<N) of the provided reference heart 

20 cycles, comparing each chosen heart cycle graph v(t;k) 
with each of the other chosen heart cycle graphs, and 
identifying a pair of graphs, v(t;kl) and v(t;k2), for 
which p 1 (kl ;k2) is a maximum value, where p 1 is a first 
numerical-valued, graph comparison metric that com- 

2 s pares two heart cycle graphs, v(t;k') and v(t;k"), and 
provides a first numerical metric value p 1 (k';k") that is a 
measure of how closely the two graphs resemble each 
other; 

(iv) forming a composite heart cycle graph, denoted v(t; 

30 kl 2), using graph values associated with the graphs v(t; 

kl) and v(t;k2); 

(v) determining when at least a first selected fraction fl 
(0.5sfl<l) of the N heart cycle graphs resemble the 
composite graph, according to a first condition pl(kl2; 

35 n)-(thrl)=Dl, with D 1 being always non-negative or 
always non-positive, where (thrl ) is a selected first met- 
ric threshold value; 

(vi) when at least the selected fraction fl of the N heart 
cycle graphs resemble each other according to the first 

40 condition, creating or providing a reference heart cycle 
graph v(t;ref)=v(t;kl2); for the reference person; 

(vii) when at least the fraction fl of the N heart cycle graphs 
do not resemble each other, incrementing the count 
index (c^c+1) to provide an incremented count index; 

45 (viii) when the incremented count index c is not greater 
than a selected maximum index value c(max), returning 
to step (iii) and repeating steps (iii)-(viii) at least once; 
and 

(ix) when the incremented count index satisfies c>c(max), 

50 interpreting this condition as indicating that it is not 
possible to create a reference heart cycle graph for the 
reference person. 

2 . The method of claim 1 , wherein said process of compar- 
ing two heart cycle graphs comprises: 

55 (iii- 1 ) forming a correlation of said heart cycle graphs, 

v(t;k=kl ) and v(t;k=k 2 ), according to a relation 


M 

corral; k2) = ^ w(l m )v(t m \ tl ) * v(L mk ; k2). 


for a selected sequence of times, t=t m (m=l, . . . , M), where 
w(t m ) is a selected non-negative weight value and * indicates 
65 a conjugation operation; 

(iii-2) computing said first metric value p 1 (kl ;k2) accord- 
ing to a relation 
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pi (kl ;k2)=corr(kl ;k2)/ {corr)kl ;kl )corr(k2;k2)} 1/2 ; and 

(iii-3) when said first metric pi satisfies a condition Ipl 
(kl;k2)[-(thrl)=Dl>0, fora selected positive threshold 
metric value (thrl), interpreting this condition as indi- 
cating that the heart cycle graphs v(t;kl) and v(t;k2) 5 
resemble each other. 

3. The method of claim 2, further comprising: 

(iii-4) whenpl satisfies lpl(kl;k2)l-(tlirl)=Dl<0,forsaid 
threshold metric value (thr), interpreting this condition 
as indicating that the heart cycle graphs v(t;kl) and to 
v(t;k2) do not resemble each other. 

4 . The method of claim 1 , wherein said process of compar- 
ing two heart cycle graphs comprises 

(iii-1) forming said first metric as a weighted sum of dif- 
ferences of said reference heart cycle graphs, v(t;k=kl ) 15 
and v(t;k=k2), according to a relation 


M 

p(kl; kl) = ^ v/(r m )|{v(r m ; *1)/ V(*1)J - (v(r„ ; *2)/ V(*2)}|'», 

m=l 

M 

v(k) = Yj lv(t»;*)l. 


for a selected sequence of times, t=t m (m=l, . . . , M), where 
w'(t m ) is a selected non-negative weight value, q is a selected 
positive number, and {v(t m ,k)} m is a selected sequence of 
heart cycle values, not necessarily associated with said can- 
didate-person; and 30 

(iii-2) where pl(kl;k2) satisfies Ipl(kl;k2)l-(thrl')s0, for 
a selected positive threshold metric value (thrl'), inter- 
preting this condition as indicating that the heart cycle 
graphs v(t;kl) and v(t;k2) resemble each other. 

5 . The method of claim 4 , further comprising: 35 

(iii-3) where pl(kl;k2) satisfies I p 1 (kl ;k2)l— (tlirl ')>0, for 

said threshold metric value (thrl'), interpreting this con- 
dition as indicating that the heart cycle graphs v(t;kl) 
and v(t;k2) do not resemble each other. 

6. The method of claim 1, further comprising providing 40 
said reference heart cycle graph by a process comprising: 

defining a maximum first metric value pl(kl,k2)asa value 
equal to max^ i .»{pl(k',k")}; and 
defining said reference heart cycle graph v(t;ref) as Avg tl 


O J O X \ ? / OKI, 

i 2 {v(t;kl), v(t;k2)}. 45 

7. The method of claim 6, further comprising choosing said 
average function Avg{v(t;kl), v(t;k2)} from the group con- 
sisting of: 

Avg{v(?;£l),v(/;A2)}={v(r;£l)+v(r;£2)}/2, (i) ^ () 

Av%{v(t;kl),v(t;k2)}=sgi\{v(t;kl)+v(t;k2)} \v(t;k\yv(t; 

&2)l 1/2 , and “ (ii) 

Avg{v(r;Xl),v(?;X2)}=v(?;A:l)+v(ci2)-sgn{v(r;Xl)+v(r; 55 

k2)Y\v(t;k\ )-v(t;k2) 1 1/2 . (iii) 

8. The method of claim 1, further comprising: 


passing said sequence of N reference graphs through a 
digital filter that (i) suppresses or removes all frequency 
contributions to said N reference graphs below a first 60 
threshold frequency f(low;l) and (ii) removes or sup- 
presses all frequency contributions to said N reference 
graphs above a second threshold frequency f(high;2) 
that is greater than the first threshold frequency f(low; 1 ). 

9. The method of claim 8, further comprising choosing said 65 
first threshold frequency f(low; 1 ) to be no greater than about 
2 Hz. 
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10 . The method of claim 8 , further comprising choosing 
said second threshold frequency f(high;2) to be no less than 
about 42 Hz. 

11 . The method of claim 1 , further comprising: 

(x) initializing a count index c' (c'(initial)=0); 

(xi) providing or receiving a sequence of N' heart electrical 
action cycles, numbered n'= 1, . . . , N', from a candidate- 
person (c-p), and forming a graph v'(t;n';c-p) represent- 
ing each heart cycle for the candidate-person; 

(xii) randomly choosing K' c-p heart cycle graphs, num- 
bered k'=l, 2, .... K' (3>K'>N')ofthe provided c-p heart 
cycle graphs, comparing each chosen c-p heart cycle 
graph v'(t;k';c-p) with each of the other chosen c-p heart 
cycle graphs, and identifying a pair of heart cycle 
graphs, v'(t;kl';c-p) and v'(t;k2';c-p), for which a value 
pl(kl';k2';c-p) is a maximum value; 

(xiii) forming a c-p composite heart cycle graph, denoted 
v'(t;kl 2';c-p), using data associated with the graphs v'(t; 
kl';c-p) and v'(t;k2';c-p); 

(xiv) determining when at least a second selected fraction 
f2 (0.5<f2<l) of the N' c-p heart cycle graphs resemble 
the c-p composite graph, according to a second condi- 
tion pl(kl2';n')-(thrT)=Dr, where (thrl') is a second 
selected first metric threshold value and Dl' is always 
non-negative or is always non-positive; 

(xv) when at least the second selected fraction £2 of the N' 
c-p heart cycle graphs resemble each other, providing a 
representative heart cycle graph v'(t;rep)=v'(t;kl2';c-p); 

(xvi) comparing the representative graph v'(t;rep)=v'(t; 
kl2';c-p) for the candidate-person with said reference 
graph v(t;ref)=v(t;kl2) for said reference person by 
comparing two graphs, v'(t;rep) and v(t;ref), using a 
second, numerical valued, graph comparison metric p2 
that compares two graphs and provides a numerical sec- 
ond metric value p2(rep;ref) that is a measure of how 
closely the two graphs resemble each other; 

(xvii) when p2(rep;ref) satisfies a third condition p2(rep; 
ref)-(thr2)=D2, where (thr2) is a selected second metric 
threshold value and D2 is always non-negative or is 
always non-positive, interpreting satisfaction of the 
third condition as indicating that the candidate-person is 
likely to be said reference person; 

(xviii) when p2(rep;ref) satisfies p2(rep;ref)-(thr2)=D2", 
where D2" is negative at least once and is positive at least 
once, interpreting this condition as indicating that the 
candidate-person is not likely to be said reference per- 
son; 

(xix) when at least the second fraction f2 of the N' heart 
cycle graphs do not resemble each other, incrementing 
the count index (c 1 — »c'+l), to provide an incremented 
count index c'(incr); 

(xx) when the incremented count index c'(incr) is not 
greater than a selected maximum index value c'(max), 
returning to step (xii) and repeating steps (xii)-(xx) at 
least once; and 

(xxi) when the incremented count index satisfies c'(incr) 
>c'(max), interpreting this condition as indicating that it 
is not possible to create a representative heart cycle 
graph for the candidate-person. 

12 . The method of claim 11 , wherein said process of com- 
paring two heart cycle graphs comprises: 

(xvi-1 ) forming a correlation of said c-p heart cycle graphs, 
v'(t;k=kl ;c-p) and v'(t;k=k2;c-p), according to a relation 
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M 

coirfl I ; k2' c — p) = ^ W kl; c — p)* v r (t m h ; k'2\ c — p). 


for a selected sequence of times, t=t m (m=l, . . . , M), where 
w'(t m ) is a selected non-negative weight value and * indicates 
a conjugation operation; 

(xvi-2) computing said second metric value p2 according 
to a relation 10 

p2(kl ;k2;c-p)=corr(kl ;k2;c-p)/{corr(kl ;kl ;c-p)corr(k2; 
k2;c-p)} 1/2 ; and 

(xvi-3) whenp2 satisfies a condition Ip2(kl;k2;c-p)-(thr2) 
=D2>0, for said second metric threshold value (thr2), J5 
interpreting this condition as indicating that said c-p 
heart cycle graphs v'(t;kl ;c-p) and v'(t;k2;c-p) resemble 
each otlier. 

13 . The method of claim 11 , further comprising: 

(xiii-4) where p2 satisfies p2(kl ;k2)-(thr2)=D2<0, for 20 
said threshold metric value (thr2), interpreting this con- 
dition as indicating that the heart cycle graphs v'(t;kl ;c- 
p) and v'(t;k2;c-p) do not resemble each other. 

14 . The method of claim 11 , wherein said process of com- 
paring two heart cycle graphs comprises 25 

(xvi-1) forming said second metric as a weighted sum of 
differences of said c-p heart cycle graphs, v'(t;k=kl ;c-p) 
and v'(t;k=k2;c-p), according to a relation 


p2(kV, k2\ c- p) = 

M' 


M' 

V(k) = Y | v'(v ; k\ c - p ) |, 40 


for a selected sequence of times, t=t m (m’ 1, . . . , M'), where 
w'(t m ,) is a selected non-negative weight value, q' is a selected 
positive number, {v'(t„,,.;k;c-p)} m is a selected sequence of 
heart cycle values, not necessarily associated with said can- 
didate-person, and * indicates a conjugation operation; and 
(xvi-2) when the second metric p2 satisfies p2(kl;k2;c-p)- 
(thr2)<0, for a selected metric threshold value (thr2), 50 
interpreting this condition as indicating that the heart 
cycle graphs v'(t;kl;c-p) and v'(t;k2;c-p) resemble each 
other. 

15 . The method of claim 14 , further comprising: 

(xvi-3) when p2 satisfies p2(kl ;k2)-(thr2)>0, for said 55 
threshold metric value (thr2), interpreting this condition 
as indicating that said heart cycle graphs v'(t;kl ;c-p) and 
v'(t;k2;c-p) do not resemble each other. 

16 . The method of claim 11 , further comprising providing 
said representative c-p heart cycle graph by a process com- 60 
prising: 

defining a maximum first metric value pl(kl,k2)asa value 
equal to max t , t „{p 1 (k',k")J ; and 
defining said representative c-p heart cycle graph v(t;rep) 
as Avg A . ltA . 2 {v(t;kl), v(t;k2)}, where Avg tl t2 {v'(t;kr;c- 65 
p), v'(t;k2';c-p)} is an average value for said heart cycle 
graphs v'(t;kl'c-p) and v'(t;k2';c-p). 


17 . The method of claim 16 , further comprising choosing 
said average function Avg{v(t;kl';c-p), v(t;k2';c-p)} from the 
group consisting of: 

Awg{v\t;k'l;c-p) ,v'(t;k2' ;c-p)}={ v' (t;k\ ';c-p)+v \t;k2 

c-p)}/2, (i) 

Avg{v , (t;kl;c-p),v'(t;k2';c-p)}=sgpi{v'(t;kl , ;c-p)+v‘(t; 

k2';c-p)}\v'(t;kl';c-pyv'(t;k2';c-p)\ u2 , and (ii) 

As/g{v\t;kV;c-p)y(t;k2';c-p)}=v\t;kV;c-p)+vXt;k2';c- 
p)-s%a{v'(t;k\ ;c-p)+v'(t;k2' ;c-p)}-\v'(t;kV ;c-pyv’ 

(t;k2‘;c-p)\ la . (iii) 

18 . The method of claim 11 , further comprising: 
estimating a heart beat rate HR, expressed in beats per 

minute (bpm) or an equivalent time rate for a heart cycle 
signature, for said reference person and for said candi- 
date-person, where the heart electrical cycle signature 
comprises a first time interval plus a second time interval 
that does not overlap the first time interval except pos- 
sibly at interval end points, where the first time interval 
ends and the second time interval for said reference 
person begins approximately at a central point, t=t(acp; 
ref) and the first time interval and the second time inter- 
val for said c-p person begins approximately at a central 
point, t=t(acp;c-p); 

providing a first normalization factor al , relating a ratio of 
a length of the first time interval for the HR for said 
reference person to a first reference length of the first 
time interval for HR=60 bpm or the equivalent time rate; 
providing a second normalization factor a2, relating a ratio 
of a length of the first time interval forthe HR for said c-p 
person to a second reference length of the first time 
interval for HR=60 bpm or the equivalent time rate; 
providing a third normalization factor (5 1 , relating a ratio of 
a length of the second time interval for the HR for said 
reference person to a third reference length of the second 
time interval for HR=60 bpm or the equivalent time rate; 
providing a second normalization factor (52, relating a ratio 
of a length of the second time interval for the HR for said 
candidate-person to a fourth reference length of the sec- 
ond time interval for HR=60 bpm or the equivalent time 
rate; 

replacing said reference graph v(t;ref) by a corrected or 
modified reference graph that is defined as 

v(f;ref)«,„. = v{f(acp);ref)-al (r(acp);ref)-Y; ref)} (Ar 
(acp);ref)). 

vh;ref) OTrr =v{r(acp);ref)+pi(r-r(acp);ref);ref)} 

(/(acp);ref)<sf); and 

replacing said representative graph v’(t;rep) by a corrected 
or modified representative graph that is defined as 

rY/jrcp)^^ v'{/(acp);t>^?)- 0 l 2 (f(acp ):c-p )-t; rep)} 

(at(acp);c-p)), 

v'(r;rep) co ^=v'{r(acp);c-/?)+p2('r-r(acp);o/?);rep)} 

(r(acp);c-p)sr). 

19 . The method of claim 18 , further comprising defining 
said first normalization factor a=al or a=a2 to be 
a={ 195. 36-0.384(HR)}/{ 195. 36-0.384(60)}. 

20 . The method of claim 18 , further comprising defining 
said first normalization factor |3=(51 or (52 to be |5={HR/60} r , 
where r is a selected fraction in a range Visr-zVi. 

21 . The method of claim 11 , further comprising: 
passing said K' c-p heart cycle graphs through a digital 

filter that (i) suppresses or removes all frequency contri- 
butions to said K' c-p heart cycle graphs below a third 
threshold frequency f(low;3) and (ii) removes or sup- 
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presses all frequency contributions to said K' c-p heart 
cycle graphs above a fourth threshold frequency f(high; 
4) that is greater than the third threshold frequency 
f(low;3). 

22 . The method of claim 11 , further comprising choosing 
said third threshold frequency f(low;3) to be no greater than 
about 2 Hz. 

23 . The method of claim 11 , further comprising choosing 
said fourth threshold frequency f(high;4) to be no less than 
about 42 Hz. 

24 . A method of authenticating or declining to authenticate 
an asserted identity of a candidate-person, the method com- 
prising: 

(i) initiating a count index c' (c'(initial)=0); 

(ii) receiving an assertion, by a candidate-person (c-p), that 
the candidate-person is identical with a specified refer- 
ence person; 

(iii) providing a sequence of N' heart electrical action 
cycles of the candidate-person, and forming a c-p heart 
cycle graph representing each provided c-p heart cycle; 

(iv) randomly choosing K' graphs, numbered k'=l, 2, , 

K' (3<K'sN) of the provided c-p heart cycles, and com- 
paring each chosen c-p heart cycle graph with each of the 
other chosen c-p heart cycle graphs; 

(v) determining whether a majority K'^ of the K' c-p heart 
cycle graphs resemble each other, according to a second 
selected numerical -valued, graph comparison metric p2 
that compares two heart cycle graphs, k'=k" and k'=k'", 
and provides a second graph comparison metric value 
p2(k";k"') for the comparison; 

(vi) when the majority K' M of the K' c-p heart cycle graphs 
resemble each other, providing a representative c-p heart 
cycle graph representing the majority K' M of the chosen 
c-p heart cycle graphs; 
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(vii) providing a reference heart cycle graph that corre- 
sponds to the identity asserted by the candidate-person; 

(viii) comparing the representative c-p heart cycle graph 
with the reference heart cycle graph, using the second 
graph comparison metric p2(rep;ref), applied to the rep- 
resentative heart cycle graph and to the reference heart 
cycle graph; 

(ix) when the value p2(rep;ref) provided by the second 

10 metric p2 satisfies p2(rep;ref)-p2(thr2)=D2, where 

p2(thr2) is a second selected threshold value and D2 is a 
real number that is always non-negative or that is always 
non-positive, interpreting this condition as indicating 
that the candidate-person is likely to be the reference 
15 person; 

(x) when the value provided by the second metric p2c does 
not satisfy the condition p2(rep;ref)-p2(thr2)=D2, with 
D2 always non-negative or D2 always non-positive, 
interpreting this condition as indicating that the candi- 

20 date-person is not likely to be the reference person; 

(xi) when the majority ¥2 M of the N' c-p heart cycle graphs 
do not resemble each other, incrementing the count 
index c' (c'— »c'+l); 

, 5 (xii) when the incremented count index c' is no greater than 

a selected maximum index value c'(max), returning to 
step (iv); and 

(xiii) when the incremented count index c' satisfies c'>c' 
(max), interpreting this condition as indicating that it is 
30 not possible to create a representative heart cycle graph 
for the candidate-person. 





